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Abstract 



We analyze the influence of spatial orientation on the optical response of 
hydrogenated silicon quantum wires. The results are relevant for the inter- 
pretation of the optical properties of light emitting porous silicon. We study 
(lll)-oriented wires and compare the present results with those previously 
obtained within the same theoretical framework for (OOl)-oriented wires [F. 
Buda, J. Kohanoff and M. Parrinello, Phys. Rev. Lett. 69, 1272, (1992)]. In 
analogy with the (OOl)-oriented wires and at variance with crystalline bulk 
silicon, we find that the (lll)-oriented wires exhibit a direct gap at k = 
whose value is largely enhanced with respect to that found in bulk silicon 
because of quantum confinement effects. The imaginary part of the dielectric 
function, for the external field polarized in the direction of the axis of the 
wires, shows features that, while being qualitatively similar to those observed 
for the (001) wires, are not present in the bulk. The main conclusion which 
emerges from the present study is that, if wires a few nanometers large are 
present in the porous material, they are optically active independently of their 
specific orientation. 

PACS number (s): 78.55 Hx, 73.20 Dx, 78.65 Gb 
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I. INTRODUCTION 



Silicon, one of the most relevant materials from the technological point of view, has 
excellent electronic properties as a semiconductor. However, in the crystalline form, it 
does not exhibit good optoelectronic properties because of the indirect gap, in the infrared 
range, between the valence and the conduction bands, and of the low probability of radiative 
electron-hole recombination. However it is possible to modify deeply such properties using 
a process of electrochemical etching which produces a material known as light emitting 
porous silicon (LEPS). This form, because of its luminescence properties, is presently the 
object of great scientific and technological attention f^. In fact, inside this material one can 
observe some wire structures: the wires are strictly interconnected with each other, and have 
a mean diameter of a few tens of Angstrom. Quantum confinement of charge carriers seems 
to be at the origin of the optical characteristics of the material, which are quite different 
from those of bulk silicon. 

In the recent past many different theoretical studies [pHTOH have been developed with 



the aim of analyzing the actual validity of the quantum confinement hypothesis in order to 
explain those new optical effects. According to the current experimental evidence, besides 
the presence of nanostructures of reduced dimensionality (wires and/or dots), the key ingre- 
dients that are necessary to induce luminescence are crystallinity and surface passivation of 
dangling bonds. In the present work we concentrate on the quantum wire model proposed 
to explain the optical properties of LEPS. We use the ab initio Molecular Dynamics (AIMD) 



method, as proposed by Car and Parrinello |pj]-|U|, to study silicon quantum wires oriented 



along the (111) crystallographic direction. The AIMD method has been successfully applied 



to a variety of different systems and problems |jT2| . Porous silicon samples grown on a silicon 



substrate in the (111) lattice direction were characterized both in the morphology and in 



the optical properties |1^-[16|. The study presented here follows a similar work developed 
for quantum wires oriented along the (100) crystallographic direction 0. This study showed 
the presence of a direct gap in the visible range and an important enhancement of the dipole 



matrix elements, which are directly related to the transition probability between the valence 
and the conduction band and, in consequence of that, to the absorption coefficient. This 
work, such as all the others quoted above, shows how the optical and electronic properties 
of silicon become deeply modified by quantum confinement. These predictions, with the 
limitations due to the Local Density Approximation (LDA) used in the AIMD framework, 
support the thesis that the origin of porous silicon photoluminescence lies in the localization 
of charge carriers inside ordered structures of reduced dimensionality. 

The present work tries mainly to ascertain if the orientation may be relevant in deter- 
mining the optical and electronic properties of the material since, a priori, it is not obvious 
that charge carriers confinement acts in a similar way along different symmetry planes. In 
section II we describe the geometry of the model used and give some computational details. 
In section III the results on the structural, vibrational and electronic properties of the model 
are discussed. The final section is devoted to the conclusions. 

II. SILICON WIRE GEOMETRY AND TECHNICAL DETAILS 

Crystalline silicon is characterized by a diamond structure, in which atoms are distributed 
in a face centered cubic (fee) lattice, each lattice site being occupied by a pair of particles 
whose chemical bond, about 2.35 A long, lies in the (111) direction. The silicon lattice 
structure along the (111) direction can be characterized by three double layers periodically 
repeated with a periodicity a-\/3, where a is the lattice parameter of bulk crystalline silicon 
(5.43 A). Fig. ^ shows the smallest wire that one can build, still preserving the crystalline 
order along the (111) direction. The unit cell, indefinitely replied in space through periodic 
boundary conditions (PBC), contains 26 silicon atoms, plus 30 hydrogen atoms which are 
necessary to saturate the dangling bonds at the surface of the wire. The geometry of the 
model is such that it exhibits a non uniform size along the longitudinal direction. In a way, 
this feature makes the model more similar to real samples of porous silicon. In fact, the 
diameter of the model has a maximum value of 7.67 A, and a mean value of 5.51 A. We 
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have also studied a larger wire whose basic structure consists of 62 silicon atoms and 42 
hydrogen atoms, with a maximum diameter of 11.72 A, and a mean value of 10.37 A. We 
note that a system consisting of 104 atoms is rather demanding on the computational side 
for ab-initio Car-Parrinello simulations. For this reason, we limited our calculations for the 
larger wire to those features which may confirm some specific properties of the thinner wire. 
In both cases the size of the supercell in the longitudinal direction corresponds exactly to 
the bulk silicon characteristic size of 9.41 A along the (111) direction. The lateral size of 
the supercell is much larger than the diameter of the wire so as to avoid interactions of the 
particles with their images. The parameters needed in the AIMD simulation were suggested 
by the previous experience on this material: we used a time step of 7 a.u., an energy cutoff 
for the plane- wave expansion of 8 Ry, an electronic mass of 300 a.u., and pseudopotentials 
of the Bachelet-Hamann-Schluter |1^ type for describing the interaction between valence 
electrons and the nuclei with core electrons. 



III. RESULTS 



A. Structural relaxation and vibrational properties 

As a first step, we wanted to verify the structural stability of the wire. To this aim, 
after the electronic ground-state relative to the initial ionic configuration was found, we let 
all the ionic degrees of freedom relax with a steepest descent technique. We notice a small 
increase of the distance of the most external atoms from the axis of the wire, of the order 
of 0.16 A, analogously to what was observed in [^. As expected, we observe no relevant 
variations of the characteristic distances, along the longitudinal direction. Then, we let the 
system evolve dynamically for a time of about 2 ps, a considerable time in AIMD, so that 
one may be confident that the system is reasonably close to thermal equilibrium. During this 
evolution the model did not exhibit any pathological behaviour: the system oscillates about 
the equilibrium configuration with no substantial atomic rearrangements. We observe that 
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the total energy (potential + kinetic) fluctuates around an average value which remains 
constant, as should be in a microcanonical AIMD simulation. The average temperature 
during this dynamical evolution is about 70 K. 

We also computed the velocity autocorrelation function and, through its Fourier trans- 
form, an indicative vibrational density of states (VDOS) of the system. We flnd a good 
qualitative agreement with the typical experimental data of the vibrational modes of Si-Si 



and Si-H bonds . From a strictly quantitative point of view, we remark that the relative 
intensity of the peaks should not be taken too rigorously. However the positions of the peaks 
are in fair agreement with the experimental data. In particular we recognize in the spectra 
of Fig. the characteristic signal of the stretching of the Si-H bond at about 2000 cm~^. 
The double peak has its origin in the presence of doubly hydrogenated silicon atoms, which 
give rise to symmetrical and antisymmetrical stretching of the bond at slightly different 
frequencies. The analysis of the oscillations of single and twin hydrogen particles shows 
that the higher-frequency peak can be ascribed to the symmetrical mode, while the lower- 
frequency signal arises from antisymmetrical stretching vibrations. We notice also, at lower 
frequencies, the typical signals of scissors mode ( 800 cmT^) of the twin hydrogen atoms 
bonded to one silicon atom, and of Si-H bending (around 600 cm~^). The weak structure 
centered at about 1150 cm~^ is likely associated with another bending vibrational mode of 
hydrogen atoms as it shows up in the motion of both single and twin passivating hydrogens. 
The low frequency noise produced by the flnite time evolution renders more uncertain the 
identiflcation of the modes below 500 crrT^ that are associated with the motion of silicon 
atoms. It is important to remark that such characteristic vibrational modes are observed. 



by infrared spectroscopy, also in real freshly prepared LEPS samples ||T8 |. 

For the larger size sample, we have also studied the relaxation towards the global energy 
minimum. It is noteworthy that in this case the mean radial relaxation, in the equilibrium 
configuration, is negligibly small if compared with that estimated for the thinner wire. This 
can be intuitively understood if one thinks that in the smaller-size model the ratio of surface 
to bulk atoms is definitely higher than in the larger wire. Such a small relaxation originates 



from the competition between the doubly hydrogenated sihcon atoms, which move shghtly 
away from the axis of the wire, and the singly hydrogenated ones, which instead move 
inward. 

B. Electronic and optical properties 

Once the minimum energy configuration has been determined, we computed the elec- 
tronic band structure of the quantum wire along the one-dimensional Brillouin zone (BZ). 
We computed the Kohn-Sham eigenvalues for six different k-points of the BZ, from the 
center (F) to the border (which would correspond to L in the bulk). The resulting band 
structure is shown in Fig. ^ for the (lll)-oriented wire with a mean diameter of 5.51 A. 
The main result is the prediction of a direct gap of 2.83 eV between the conduction and the 
valence band at the zone center. It is also interesting to notice that the second conduction 
state still shows a signature of the band structure of bulk silicon, i.e. a minimum located 
at about 0.2. As to the gap, we recall that Buda and co-workers found, for a (100) -oriented 
quantum wire with a diameter of 7.60 A, a value of 2.33 eV [Q. The difference is to be 
ascribed, as we shall see later, to the different mean diameter of the two models, more than 
to the different symmetry of the wires. We recall that in bulk silicon the minimum of the 
conduction band along the F-L direction falls just at L, with a gap of 2.15 eV. 

For the larger wire we find an energy gap of 1.90 eV at the center of BZ which, consistently 
with the quantum confinement hypothesis, is larger than that in bulk silicon, but smaller 
than the value for the thinner wire. The values of the two energy gaps, together with those 
of the (100) models and of bulk silicon as obtained in LDA, are plotted in Fig. |^. It is of 
striking evidence the linear behaviour of the energy gap as a function of the inverse size of 
the models. It is remarkable that both the (111) and (100) computed values lie on the same 
straight line, suggesting the independence of such property on the specific orientation of 
the wires. In the same figure we also show the theoretical prediction coming from effective 
mass arguments {Eg ~ ^), which should hold close to the bulk limit. It appears that the 
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matching between these two different behaviours should occur in the region between 2 and 
5 nm. 

The presence of a direct gap, and thus the possibihty of zero-phonon transitions, is 
not a sufficient condition for ensuring good photoluminescence properties. The probability 
of an electronic transition as a result of the absorption of a photon is given by the Fermi 
golden rule, which in the dipole approximation is | (V'ckle- pj'i/'ijk) P where -i/^ck and '0t,k are the 
conduction and valence states involved in the electronic transition, while e is the polarization 
vector and p the momentum operator. We determined all the matrix elements between 
Kohn-Sham states of the conduction and valence bands whose difference in energy falls 
in the visible range, or slightly above. Table | and Table || give the values of the most 
relevant of these dipole-matrix elements for the smaller and the larger (lll)-oriented wires. 
We note the strong difference between the values relative to the longitudinal and transverse 
polarizations. Moreover, the quantum confinement of the electrons mixes up the states, thus 
allowing for transitions that are forbidden in the bulk by selection rules. The dipole-matrix 
element between the top of the valence band and the bottom of the conduction band is not 
zero because of the mixing, even though not very large. 

We computed the imaginary part of the dielectric function, that is related to the absorp- 
tion coefficient, through the formula 

= E / ^ K^e.Je ■ - E^iK) - M (1) 

m"=u;"= V 7^ Jbz 2tc 

where and V are the length of the wire and the volume of the supercell, respectively. The 
results for the smaller model are shown in Fig. |^. In the construction of this quantity we 
have considered only the dipole-matrix elements corresponding to energy differencies up to 4 
eV. We observe that the component related to an external field parallel to the wire direction 
is much stronger than the perpendicular one, and this difference becomes larger and larger 
going to lower energies. Probably, the most interesting result is the presence, in the parallel 
component, of a small signal at the band gap (2.83 eV) which is completely absent in the 
perpendicular one. This is a clear effect of the mixing of the bulk states induced by the 
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reduced dimensionality. This result seems to be peculiar of the (111) wire, since in the (100) 
models the dipole-matrix element between the top of the valence band and the bottom 
of the conduction band was negligible in both the components. 

The plot of the dielectric function for the larger-sized wire, that is shown in Fig. |^, 
exhibits also a different behaviour in the parallel and the perpendicular components. We 
notice also in this case a clear difference in intensity between the parallel and the perpendic- 
ular components at higher energies. The most remarkable feature is the presence, below ~ 
2.7 eV, of three peaks in the longitudinal component that are not obseved in the transverse 
one. In particular we observe again a significative signal in correspondence of the energy 
gap (1.90 eV), whose intensity is lowered with respect to the smaller wire. This is consistent 
with the fact that in the bulk limit this matrix element should go to zero. 

IV. CONCLUSIONS 

We have presented an electronic structure calculation, within the LDA, of two (111)- 
oriented silicon quantum wires of different size. The comparison of the present results with 
those obtained for the (OOl)-oriented models shows a substantial analogy of the electronic 
and optical properties of silicon quantum wires oriented along different crystallographic di- 
rections. In particular we find a direct energy gap at the zone center whose magnitude 
decreases for larger wires consistently with the quantum confinement hypothesis. Quite 
surprisingly we find that the values of the gap for all the (111) and (001) wires studied 
are excellently fitted as a linear function of the inverse diameter. Therefore this quantity 
seems to be independent on the orientation and geometrical details, and is affected only by 
the width of the wire. Again, the present calculations support the idea that the photolu- 
minescence properties of porous silicon can be ascribed to charge carriers confinement in 
nanostructures. We find, similarly to the (001) wires, a strong anisotropy in the imaginary 
part of the dielectric function. These results indicate that, if the LEPS structure consists of a 
network of randomly oriented wires , all of them will contribute to the photoluminescence 
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independently on their specific orientation. 
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FIGURES 

FIG. 1. Silicon quantum wire oriented in the (111) direction: side view (panel a) and top 

view (panel b). The large balls represent silicon atoms while the small balls indicate the hydrogen 
atoms passivating the dangling bonds at the surface of the wire. 

FIG. 2. Vibrational spectrum for a (lll)-oriented wire with an average diameter of 5.51 A. 
We rescaled of a factor 5 the feature at about 2000 cm~^ in order to emphasize the double peak 
(see text). 

FIG. 3. Electronic band structure for the smaller (lll)-oriented wire along the one-dimensional 
Brillouin zone. We express k in unit of where c is the periodicity along the direction of the 
wire. Only few bands close to the energy gap are plotted. 

FIG. 4. Computed energy gap as a function of the inverse wire diameter. We show both the 
(111) wires (diamonds) and (100) wires (crosses) computed gap. The full line is the linear fit to this 
points. The dashed curve represents the effective mass theory prediction. The LDA bulk silicon 
gap is represented by the square on the vertical axis. 

FIG. 5. Imaginary part of the dielectric function e2(E) for the thinner wire. The plot is on a 
semilogarithmic scale. The full line represents ef, while the dashed line refers to 62 ■ 

FIG. 6. Imaginary part of the dielectric function e2(E) for the larger wire. The plot is on a 
semilogarithmic scale. The full line represents e|, while the dashed line refers to 62 ■ 
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TABLES 

TABLE I. Square modulus of the transverse and longitudinal dipole-matrix elements for the 

smaller wire in unit of Bohr~'^. The top of the valence band and the bottom of the conduction band 
are labelled by indices 67 and 68 respectively. Ec — is the difference between the conduction 
and valence Kohn-Sham eigenvalues. 
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TABLE IL Square modulus of the transverse and longitudinal dipole-matrix elements for the 
larger wire in unit of Bohr~^. The top of the valence band and the bottom of the conduction band 
are labelled by indices 145 and 146 respectively. E^ — E^ is the difference between the conduction 
and valence Kohn-Sham eigenvalues. 
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